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ABSTRACT 


Exact  dipole  distributions  in  closed  form  are  not  known 
for  any  explicitly-defined  bodies  in  uniform  translation  under  a  free 
surface.  Exact  dipole  distributions  are  known  for  some  special 
families  of  bodies  in  uniform  translation  in  an  infinite  fluid.  One 
of  these  families  consists  of  the  lenticular  cylinders,  which  are  of 
interest  in  connection  with  ships  because  they  can  have  any  preassigned 
entrance  angle.  The  exact  and  linearized  dipole  distributions  generat¬ 
ing  lenticular  cylinders  in  an  infinite  fluid  are  compared  for  various 
values  of  a  fullness  parameter.  These  dipole  distributions  are 
approximations  to  the  exact  dipole  distribution  generating  a  (vertical) 
lenticular  cylinder  in  the  presence  of  a  free  surface.  Each  of  these 
dipole  distributions  has  a  wave  resistance,  the  wave  resistance  of  the 
linearized  dipole  distribution  being  Michell's  wave  resistance  for  the 
lenticular  cylinder.  It  is  found  that  for  practical  values  of  a  full¬ 
ness  parameter  the  two  approximate  dipole  distributions  are  not  very 
close,  and  for  Froude  numbers  around.  .4  the  wave  resistances  may 
differ  by  as  much  as  40%.  This  difference  is  larger  than  the  uncertainty 
in  the  "measured"  wave  resistance,  and  it  is  proposed  that  experiments 
be  conducted  to  determine  the  relative  merits  of  the  two  approximate 
dipole  distributions  in  predicting  wave  resistance.  If  Michell's  dipole 
distribution  for  an  arbitrary  form  is  modified  in  an  ad  hoc  but  natural 
way  suggested  by  the  results  for  lenticular  cylinders,  a  comparison 
with  existing  measured  wave  resistance  can  be  carried  out.  Similar 
remarks  apply  to  the  asymptotic  behaviour  of  the  wave  resistance  at 
low  speed. 
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TABLE  OF  NOMENCLATURE 

Aj  cross  sectional  area  of  the  cylinder 


1/2 

L/2 

A  =  ^  ( 

%  c  ) 

a(x)dx  =  -  j  a(x)dx 

-1/2 

-L/2 

a,  one  half  the  width  of  the  cylinder  in  the  x-direction 

(see  Figure  2) 

C  ,  wave  resistance  coefficient 

w’ 

EIF 

C  ,  wave  resistance  coefficient  of  the  dipole  distribution 

w 

which  represents  the  body  exactly  in  an  infinite  fluid 

(and  related  to  f^(x)) 

M 

C  ,  wave  resistance  coefficient  of  Michell's  dipole  distri- 

W 

bution  (and  related  to  g  (x)) 

€ 

c,  cylinder  velocity 

f,  Froude  number  =  c/  / gL 
F  =  1/f^ 

F^,  coefficient  of  (l-4x^)€^  in  a  power  series  expansion 

of  f^(x) ,  see  Equation  (10) 

t.(x)  =  -  ^  OeIf'") 

G,  Green's  function 

G^,  coefficient  of  (l-4x  ) e"  in  a  power  series  expansion  of 

g^(x) ,  see  Equation  (9) 

g,  gravitational  acceleration 

2 

~  ~  cross  sectional  shape  of  the  lenticular 

cylinder 
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h(x)  , 


a  function  proportional  to  a  dipole  distribution,  such 

1/2 

that  j  h(x)dx  =  1. 

-1/2 


radius  of  the  circle  forming  the  cross  section  of  the 
lenticular  cylinder 


L  =  2a,  see  Figure  1, 


an  integer  greater  than  zero  such  that  mp  ^  1 
parameter  defining  the  shape  of  the  cylinder,  see 
Figure  2, 
wave  resistance 


wave  resistance  of  a, 


/^  A  A 


wave  resistance  of 

local  radial  coordinate,  see  Figure  8 
fluid  velocity 

1 

complex  velocity  potential 

normalized  rectangular  coordinates,  see  Figure  1  and 
Appendix  B. 

unnormalized  rectangular  coordinates 

the  x-coordinate  of  the  m^^  singularity  in 

Bessel  function  of  the  second  kind,  order  zero 

the  z  coordinate  of  the  m^^  singularity  in  ^ 

local  angular  coordinate,  see  Figure  8 

average  half-beam  of  the  shape,  divided  by  the  length, 

i.e. ,  c  =  A/2L^ 


1/2 

Cjr  =  /  f^(x)dx 

-1/2 
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1/2 


g^(x) dx 


C 


h 


e, 


e  +  iTi 

bipolar  coordinate  defined  by  Equations  (1)  and  (2) 
strength  of  the  discrete  dipole  source 
bipolar  coordinate  defined  by  Equations  (1)  and  (2) 


e 


o 


р. 

A 

с,  a, 


'^EIF’ 


"M’ 


density 

dipole  distribution;  ct(x)  =  ct(x) 

dipole  distribution  which  is  exact  in  an  infinite  fluid 
Michell's  dipole  distribution 
velocity  potential,  see  Appendix  B 

approximate  ^  in  the  vicinity  of  the  discrete  dipole 
source . 
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1 .  INTRODUCTION 

In  1898  J.H.  Michell  published  a  paper  entitled  ’’The  Wave 
Resistance  of  a  Ship",  in  which  he  gave  a  formula  for  the  approximate 
wave  resistance  of  a  thin  ship^^^.  It  was  found  later^^^’  by  making 

a  formal  expansion  of  the  complete  nonlinear  problem  in  terms  of  a 
small  thickness  parameter  such  as  the  beam,  that  Michell' s  formula 
correctly  predicts  the  leading  term  in  the  wave  resistance.  The  formal 
expansion  procedure  can  be  used  to  find  higher-order  approximations  to 
the  solution  of  the  exact  nonlinear  problem  but  the  calculations  are 
complex  and  have  never  to  our  knowledge  been  carried  out  for  any  ship 
form. 

Michell 's  result  is  obtained  by  linearizing  the  boundary 
conditions  on  the  ship  and  on  the  free  surface.  If  we  satisfy  one  of 
these  boundary  conditions  to  a  higher  order  than  Michell  did  we  will 
obtain  a  new  approximate  solution  to  the  exact  nonlinear  problem.  This 
new  solution  will  agree  with  Michell  (and  hence  be  exact)  to  the  first 
nonvanishing  order  in  the  thickness  parameter.  The  second  terms  of 
Michell 's  solution  and  the  new  solution  will  disagree  and  both  are 
presumably  incorrect.  ^'Jhen  applied  to  any  particular  case  one  of  these 
solutions  will  be  numerically  more  accurate  than  the  other.  It  may  be 
that  our  solution  is  more  accurate,  in  which  case  it  would  be  useful. 

We  cannot  prove  that  it  is,  but  we  do  show  that  our  solution  has 
desirable  properties. 

The  problem  cf  satisfying  the  boundary  condition  on  the  ship 
to  higher  order  arises  also  in  problems  without  free  surface.  In  fact 
the  uniform  flow  past  a  lenticular  cylinder  in  an  infinite  fluid  is 
kno\^m  exactly.  This  exact  solution  and  the  corresponding  linearized 
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1 ■  INTRODUCTION 

In  1898  J.K.  Michell  published  a  paper  entitled  "The  Wave 
Resistance  of  a  Ship",  in  which  he  gave  a  formula  for  the  approximate 
wave  resistance  of  a  thin  ship^^^.  It  was  found  later^^^’^^^  by  making 
a  formal  expansion  of  the  complete  nonlinear  problem  in  terms  of  a 
small  thickness  parameter  such  as  the  beam,  that  Michell' s  formula 
correctly  predicts  the  leading  term  in  the  wave  resistance.  The  formal 
expansion  procedure  can  be  used  to  find  higher-order  approximations  to 
the  solution  of  the  exact  nonlinear  problem  but  the  calculations  are 
complex  and  have  never  to  our  knowledge  been  carried  out  for  any  ship 
form. 

Michell 's  result  is  obtained  by  linearizing  the  boundary 
conditions  on  the  ship  and  on  the  free  surface.  If  we  satisfy  one  of 
these  boundary  conditions  to  a  higher  order  than  Michell  did  we  will 
obtain  a  new  approximate  solution  to  the  exact  nonlinear  problem.  This 
new  solution  will  agree  with  Michell  (and  hence  be  exact)  to  the  first 
nonvanishing  order  in  the  thickness  parameter.  The  s:!cond  terms  of 
Michell' s  solution  and  the  new  solution  will  disagree  and  both  are 
presumably  incorrect.  When  applied  to  any  particular  case  one  of  these 
solutions  will  be  numerically  more  accurate  than  the  other.  It  may  be 
that  our  solution  is  more  accurate,  in  which  case  it  would  be  useful. 

We  cannot  prove  that  it  is,  but  we  do  show  that  our  solution  has 
desirable  properties. 

The  problem  of  satisfying  the  boundary  condition  on  the  ship 
to  higher  order  arises  also  in  problems  without  free  surface.  In  fact 
t'le  uniform  flow  past  a  lenticular  cylinder  in  an  infinite  fluid  is 
knn\,m  exactly.  This  exact  solution  and  the  corresponding  linearized 
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solution  are  discussed  in  Section  2,  particular  attention  being  paid  to 
the  exact  and  linearized  dipole  distributions  generating  these  flows. 

In  Section  3  we  use  the  results  of  Section  2  to  construct  a 
new  approximate  solution  to  the  problem  of  a  vertical  lenticular  strut 
in  a  uniform  stream  with  free  surface.  This  solution  has  the  following 
properties  t 

1)  It  agrees  with  Michell  and  is  therefore  correct  to  the 
first  nonvanishing  order  in  the  thickness  parameter 
(for  any  fixed  Froude  number). 

2)  It  is  correct  to  the  first  nonvanishing  order  as  the  Froude 
number  tends  to  zero  (for  any  fixed  thickness  parameter), 

Michell *s  solution  does  not  have  this  property, 

3)  It  is  correct  to  the  first  nonvanishing  order  as  the  depth 
of  the  observation  point  goes  to  infinity  (for  any  fixed 
Froude  number  and  thickness  parameter)  Michell* s  solution 
does  not  have  this  property. 

Note  that  while  the  new  solution  has  some  properties  not  enjoyed 
by  Michell 's  it  is  not  claimed  to  be  correct  to  higher  order  in  the  thick¬ 
ness  parameter 

In  Section  4  we  compute  a  new  wave  resistance  for  vertical 
lenticL'lar  struts  using  the  new  approximate  solution  of  the  flow  problem. 
The  values  of  wave  resistance  obtained  by  the  two  methods  differ  by  as 
much  as  407c  for  realistic  values  of  the  thickness  parameter  and  for 
Froude  numbers  around  .4  This  difference  is  larger  than  the  uncertainty 
in  experimentally  determined  wave  resistance,  and  it  is  proposed  that 
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experiments  be  carried  out  to  determine  the  relative  merits  of  the  two 
approximations  for  predicting  wave  resistance.  It  is  also  shown  that 
the  new  approximation  leads  to  a  different  low- speed  asymptotic  behaviour 
for  the  wave  resistance. 
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2.  SHAPE  vs,  DIPOLE  DISTRIBUTION  FOR  A  LENTICULAR  CYLINDER  IN  AN  INFINITE 
FLUID 

In  an  infinite  fluid  the  cross-sectional  shape  g^(x)  of  the 
infinite  cylinder  gives  the  linearized  dipole*  distribution  -  g  (x) . 

^TT  € 

The  fullness  parameter  e  is  the  average  half-beam  of  the  form  divided  by 
its  length.  The  exact  dipole  distribution  generating  this  shape  is  defined 
to  be  -  f^(x)  o  For  thin  bodies  we  expect  g^Cx) -^f ^(x)  ,  in  accordance 

with  linearized  theory.  Since  the  cylinder  and  the  fluid  are  assumed  to 
be  infinite,  the  dipole  distribution  is  independent  of  the  y-coordinate, 
see  Figure  1,  and  we  can  confine  our  discussion  to  the  two-dimensional 
problem  in  the  xz-plane,  see  Figures  2A,  2B. 

Let  us  introduce  the  bipolar  coordinates  t); 


(1) 


A  ^  a  sinh  ri 
^  cosh  T]  -  cos^ 


(2) 


A 

Z  = 


a  sin  i _ 

cosh  T]  '  cos^ 


The  shape  of  the  lenticular  cylinder  (formed  by  the  intersection 
of  two  circular  cylinders)  is  given  by  ^  .  The  complex  potential 

w  for  such  a  cylinder  moving  in  the  negative  x-direction  with  velocity  c 
(or  for  a  stationary  cylinder  in  a  uniform  stream  moving  in  the  positive 
X'direction  with  velocity  c)  is  given  by^^^ 


(3)  w  =  ica(|  cot  ^  -  cot  ■^) 

where  '  =  ^  +  It]  The  velocity  potential  is 


*  Dipoles  in  an  infinite  fluid,  without  free  surface. 

Relationships  between  normalized  and  unnormalized  quantities  are  given 
in  Appendix  B . 
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t  -  ca(- 


sinh 


cosh  -  cos  — = 
P  P 


2a -L. 


The  parameter  p  is  related  to  €  by 


r  «  Tr(2-p)  •{  sin  p-rr  „  A 
Vd;  e  2  TTI 

8  sin  ^  2'- 

where  A  is  the  cross  sectional  area  of  the  cylinder.  Figure  3  shows  e 
vs.  p  for  1  ^  P  ^  2.  For  a  circular  cylinder  p  “  1  and  e  “  ^  ;  for  a 
strip  p  =  2  and  e  =  0.  For  p  <  1  the  shape  of  the  cylinder  is  shown  in 
Figure  2B. 

If  0  ^  p  ^  1  we  will  call  the  cylinder  a  supercircular  cylinder, 
and  if  1  <  p  ^  2  we  will  call  it  a  subcircular  cylinder.  For  a  subcircular 
cylinder  the  dipole  distribution  is  confined  to  the  x-axis  (or  the  xy-plane 
in  three  dimensions),  while  for  p  <  1  there  are  additional  discrete 
dipoles  on  the  z-axis  (i.e,,  lines  of  dipoles  passing  through  the  z-axis 
parallel  to  the  y-axis) .  In  the  rest  of  this  section  we  discuss  only  sub- 
circular  cylinders,  1  <  p  ^  2;  supercircular  cylinders  are  discussed  in 


Appendix  C . 


The  function  f^(x)  is  related  to  f)  by 


df^(x) 


I  z“0+ 


>  I  ^1-^7  • 


Using  Eq.  (4)  and  the  relationship  between  normalized  and  unnormalized 
potentials  (see  Appendix  B)  we  find 


f^(x) 


-  -  sin  (l-4x^)^/P 


(l+2x)^/P+(l-2x)^^P-2(l-4x  )^/Pcos  ^ 
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li(;  is  ^'iven  by 


(8)  gg(x)  ~  (  \/csc^  ^  -  4x^  +  cot 


For  small  €  we  can  express  f^(x)  and  g^(x)  as  a  power  series 


in  €, 


(9)  g.(x)  =  (l-'-x^)  51  G  (x)€’^ 

_ -ill 


n=l 


(10)  f  (x)  =  (l-4x^)  Z  F^(x)e’^ 

t  _ 1  II 


n*l 


where  from  Equations  (5),  (7),  and  (8),  we  find 


Gi  -  Fi  = 


G 


2 


0 


T-,  _  18  r  1  .  1  /l-2xx  1 

F2  -  V 

etc 

Thus  if  e  -►  0  (p  ->  2)  we  see  that  [g^(x)  -  f^(x)  ]  ->  0  for 
fixed  X,  since  G^  =  Fj^.  The  functions  g^(x)  and  f^(x)  are  shown  in 
Figures  4-7  for  c  =  .025,  .05,  .1,  and  ,3  .  The  fractional  difference 
of  the  two  functions  is,  for  fixed  x  and  small  e  , 


g^(x) 


f,(x) 


12 

TT 


[1  -f 


e  + 


*  The  validity  of  such  an  expansion  is  examined  in  Appendix  A, 
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v.nn’ ch  hccomoc  Infj.nj  tc  for  x  -  1  1/2,  no  matter  how  small  e  may  be. 
bcl  pvioin;  oi!  the  I'anctloris  at  x  =  1'  1/2  is  also  quite  different: 

(IT)  f  (x)''' - sin  ^  .  (7r6)^''^P  as  6  =  i  -  x  ->  0 

(12)  -  “  ban  ^  •  [ttS]  . 

T’^e  tlifrercncc  is  apparent  in  Fip,ures  4-7. 


The 
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3»  A  NEW  APPROXIMATE  FREE- SURFACE  FLOW  PAST  A  VERTICAL  LENTICULAR  CYLINDER 


Suppose  a  lenticular  cylinder  is  placed  upright,  as  shown  in 
Figure  1,  in  a  uniform  flow  with  linearized  free  surface.  The  dipole 
distribution  generating  the  resulting  mathematically  well-defined  flow 
is  unfortunately  not  known.  The  approximate  flow  proposed  by  Michell^^^ 
and  others  is  the  one  generated  by  the  dipole  distribution** 


aj^(x)  = 


^  g^(x)  for  lx  I  ^  I  ,  y  ^  0 


where 


z(x,y)  =  t  g^(x) 


is  the  equation  of  the  lenticular  cylinder.  That  this  is  the  first  non¬ 
vanishing  term  in  an  expansion  of  the  exact  nonlinear  solution  in  powers 

[2  1 

of  €  was  assumed  by  Michell  and  verified  formally  by  Stoker  and  Peters  ■' 


and  Wehausen 


distribut ion 


We  propose  the  approximate  flow  generated  by  the  dipole 


(14) 


=  '  f^(x)  for  |xl  ^  Y  >  y  ^ 


*  We  have  adopted  the  convention  v  =  +  grad  ^  ,  so  that 

222 

G(x, y , z ; x' , y ’ , z ‘ )= [ (x-x* )  +(y-y’)  +(z~z')  ]  +  regular  function 

is  the  velocity  potential  at  x,y,z  of  a  unit  sink  at  x',y’,z’  in  a  uniform 
stream  under  a  free  surface,  and  -  is  the  potential  of  a  unit 

x-directed  dipole 


**  Dipoles  sat  is f ying  die  linearized  free  surface  condition 
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To  the  first  nonvanishing  term  in  c  we  have 


(15) 


as  we  see  from  Eq„  (13),  (14)  and  Eq.  (9),  (10)  et  seq.  Also,  ffgjp(x) 
rs  exact  in  an  infinite  fluid,  as  suggested  by  the  subscript  EIF,  in 
the  sense  that  asymptotically  for  large  y  it  generates  the  exact  two- 
dimensional  flow  past  the  lenticular  cylinder.  Asymptotically  for  large 
y  0^  generates  the  linearized  two-dimensional  flow  past  the  lenticular 
cylinder. 

Finally,  if  the  Froude  number  f  tends  to  zero,  for  any  fixed 
value  of  the  thickness  parameter  e>  the  flow  defined  by  approaches, 

for  y  2  exact  two-dimensional  flow  past  the  lenticular  cylinder, 

Che  positive  images  required  to  satisfy  the  limiting  boundary  condition 

-  0  being  supplied  by  the  free- surface  dipole  potential.  Thus  the 

oy 

new  approximate  solution  to  the  mathematically  well-defined  problem  of 
a  lenticular  strut  in  a  uniform  flow  with  linearized  free  surface 
condition  has  advantages  net  enjoyed  by  Michelles  solution  to  the  same 
problem  The  suitability  of  this  model  for  analyzing  full  forms  may  be 
questioned,  since  for  finite  Froude  number  full  forms  are  inconsistent 
with  small  waves  On  the  other  hand  the  new  solution  becomes  exact  for 
low  Froude  numbers  Asymptotics  aside,  the  practical  value  of  the  new 
solution  will  be  determined  by  the  accuracy  of  its  numerical  predictions 
for  practical  forms  operating  at  practical  Froude  numbers 


% 
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^ _ A  NEW  APPROXIMATE  WAVE  RESISTANCE  FOR  A  VERTICAL  LENTICULAR  CYLINDER 

Suppose  we  have  a  dipole  distribution  a(x)  on  the  region 
1^!  y  >  0  Within  the  mathematically  well-defined  theory  of  flows 

catisfying  the  linearized  free  surface  condition  the  wave  resistance  of 

[  5  ] 

thl:  dipole  dlscribution  is  given  exactly  by 

1/2  1/2 

{  R  =  Sr^pF^  (  (  c(x)c'(x')y  (F|x'X*l  )dx  dx* 

Vi  ;'  /  •  O 

-1/2  -1/2 

where 

•:«(x)  c<xL)j 

p  is  the  density  of  the  fluid  , 

r  =  ^2  ^  f  -  1/  =  Froude  number  ^ 

c 

1,  IS  the  hessel  function  of  the  second  kind,  order  zero  . 

c  ^h 

i  f  we  ler 

A  h  -  r  e 

1/2  1/2 

^  hvx)  dx=i  and  f  c(x)dx 

W'2  1/2 

we  Car.  introduce  a  wave  resistance  coefficient 


1/2 

1/2 

(  t '  1 

c  = 

w 

2  ^ 

1 

* 

j  h ( x)h vx M Y , (F 1 X -  X M  ) dx  dx 

1/2 

-1/2 

In  term:  of  which  the  wave  resistance  becomes 

.;i8'  R  =  4  pc  rc 

w  2  ^  h  w 
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where 


If  g^(x)  defines  the  shape  of  a  subcircular  lenticular  cylinder, 
then  Michell‘s  approximate  dipole  distribution  is 

(19)  o-j^(x)  =  -  ^  8-(^)  1^1  ^  7  »  y  ^  0.  (see  Equation  (13)), 

A  M 

and  Micheli's  wave  resistance  for  the  lenticular  cylinder  is  obtained 
by  substituting  crj^(x)  for  a(x)  in  Equation  (16)  „  For  comparison  we 
propose  the  dipole  distribution 

(2Q)  ^  Ul  ^  7  .  y  2  0* 

avif 

and  the  corresponding  wave  resistance  obtained  by  substituting 

'’^If(x)  for  o(x)  in  Equation  (16).  The  properties  of  the  dipole  distri- 
Dution  CpFp  and  the  corresponding  flow  were  discussed  in  Section  3. 

We  n^w  compare  with 

OpFF  associated  with  a  given  lenticular  cylinder  will  have  different 

values  of  (  or  Cj^)  ,  we  cannot  compare  the  values  of  but  we  mtist 

2 

c  K|,,,re  values  of  (see  Equation  (18)),  which  is  done  in 

Figures  8-10 

For  Michell's  dipole  distribution  we  have 

EIF  means  Exact  in  an  Infinite  Fluid 


^EIF 
%  • 


Since  the  two  dipole  distributions  and 
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y" 

J  g^(x)dx, 

-1/2 


s 


while  for  CT^jp(x)  we  have 


^h  =  L 


i  J 

-1/2 


'  Li7 


2  M  2  EIF 

The  eraphs  show  e  C  and  e^C  vs.  the  Froude  number  f  for  cylinders 

g  w  f  w 

with  p  =  1.905,  1.812,  and  1.635.  The  e  defined  in  Equation  (5)  is  e  , 

O 


2 

Since  2L  e  is  the  cross-sectional  area  of  the  lenticular  cylinder,  while 

2 

2L  €£  is  the  dipole  moment  of  the  distribution  the  relation 

6  9 

=  -  —  €^+-‘-  indicates  that  for  small  e  the  dipole  moment  exceeds 

,3  f  tt  f  8 

the  volume.  This  is  also  true  for  a  circular  cylinder  [eg=7r/8,  e£=Tr/4]. 
Tlicse  results  are  examples  of  a  theorem  of  Landweber  and  Yih^^^,  which 
states  that  added  mass  +  buoyancy  =  (aipole  moment];  since  the 

buoyancy  is  the  density  times  the  volume,  and  the  added  mass  is  positive, 
Lhe  dipole  momenL  always  exceeds  the  volume.  It  follows  that  the 
volumetric  portion  of  the  discrepancy,  contained  in  and  e^,  is  alv/ays 
in  the  same  direction.  Michcll's  theory  tends  to  underestimate  the  wave 
resistance  because  it  underestimates  the  volume.  O/er  most  of  the  Froude 
number  range  in  Figures  8-10  the  portion  of  the  discrepancy  contained  in 
th.c  C  '  s  is  n  .uc  same  direction. 

\'7 

Figures  8-10  reveal  a  considerable  difference  in  the  values  of 
umve  resistance  computed  by  the  two  methods.  This  difference  is  discussed 
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A  pT]? 

further  below.  We  close  this  section  by  remarking  that  has  the 

fnllo-wing  property  V7hieh  follows  from  statements  made  in  Section  3:  To 
the  first  nenvanishing  term  in  e  we  have,  for  fixed  Froude  number, 


(22)  rJJ  =  , 

and  both  are  correct  to  this  order.  For  fixed  e  we  would  like  to  conclude 
that  to  the  first  nonvanishing  term  in  the  Froude  number  f,  as  f  0, 

A  p  jp  ^  W, 

is  correct  while  R^  is  not,  since  cr^jp  and  have  the  analogous 
properties.  We  have  not  been  able  to  supply  a  derivation,  since  certain 
apparently  necessary  information  about  the  manner  in  which  cJpjp  becomes 
exact,  as  f  -»>  0,  is  not  known  to  us.  However,  we  can  compare  the 
asymptotic  expansions  of  R.^  and  ^EIF 

as  f  ♦  0,  namely 


(23) 

(24) 


These  asymptotic  expansions  are  found  by  methods  used  in ^  ^  ^  and  described 

f  3  ] 

also  in  .  We  can  see  that  Equations (23)  and  (24)  are  different,  but 
decreasingly  so  as  €  -»■  0  (p  -►  2)  .  In  fact  if  we  let  p  =  2  (l-r))  ,  q  0, 
and  expand,  Equat ions (23)  and  (24)  both  yield 


.o.x  1  2.^  T?-  2  8f^  4  .  v2 

(23)  ^  '  ’  ~ir  3 

There  is  nothing  to  be  gained  by  comparing  Equations  (23)  and  (24)  since 

both  are  inadequate  for  practical  Froude  numbers.  However,  one  can  find  the 

oscillatory  terms  which  come  next  in  the  asymptotic  expansions,  at  which 

point  a  comparison  will  provide  valuable  information. 
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5.  CONCLUSIONS 

9  9 

Figures  4,  5,  6,  and  7  compare  f  =  -  g  “  ■  —  Ow 

for  a  lenticular  cylinder.  Figures  4  and  5  are  the  most  important  from 
a  practical  point  of  view  since  practical  values  of  e  (based  on  waterplane 
section)  are  in  the  range  .025-. 050  .  The  discrepancy  in  this  range  can 
be  characterized  as  follows: 

a)  In  the  central  80%  of  the  form  f^Cx)  is  larger  than  g^(x), 
the  discrepancy  being  maximum  at  the  center.  The  integrals 
are  related  as  in  Equation  (21) . 

b)  Near  the  ends  fg(x)  vanishes  to  a  different  order  than 
g^(x),  as  shown  in  Equations  (11)  and  (12). 

Figures  8,  9,  and  10  compare  the  wave  resistance  of  a^jp(x)  and 

aj^(x)  .  Figures  8  and  9  are  the  most  interesting  from  a  practical  point 
of  view.  The  numbers  were  computed  using  a  digital  computer  program 
whose  accuracy  decreases  with  Froude  number,  and  for  that  reason  no 
results  are  presented  for  f  <  .3  .  The  discrepancy  in  wave  resistance  is 
not  small  compared  to  typical  discrepancies  between  calculated  and  measured 
wave  resistance,  and  is  large  compared  to  the  uncertainty  in  measured  wave 
resistance.  The  discrepancy  appears  to  change  sign  as  a  function  of 
Froude  number.  It  would  be  interesting  to  test  lenticular  ship  models  and 
compare  measured  wave  resistance  with  the  wave  resistances  of  Oj^(x)  and 
OeifC’^)  •  It  would  be  significant  and  useful  if  a  consistently  better 
agreement  were  found  using  <^£jp(* **)  • 


*  Typical  values  of  c  are  .04  for  a  destroyer,  .048  for  a  tanker,  based 
on  the  waterplane  section. 

**  All  modified  for  finite  draft. 
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The  results  derived  above  for  a  special  class  of  forms  can  be 
used  to  construct  new  dipole  distributions  for  general  forms  having  a 
finite  entrance  angle.  This  can  be  done  at  each  depth  in  an  ad  hoc 
fashion  by  multiplying  Michell's  dipole  distribution  by  f^^^^ (x)/g^^y^ (x) , 
v/here  e(y)  defines  the  lenticular  cylinder  having  the  same  entrance  angle 
as  the  given  form  at  the  depth  y.  This  modification,  if  applied  to  forms 
for  which  model  resistance  data  are  already  available,  would  make  possible 
the  evaluation  of  a  potential  improvement  in  wave  resistance  calculation 
procedures  without  performing  additional  model  tests. 

and  predict  different  asymptotic  behaviour  of  the  wave 

resistance  at  low  speed.  It  would  be  interesting  to  conduct  model  tests 
on  lenticular  struts  to  evaluate  the  relative  merits  of  the  two  theories. 
Hov/ever,  if  we  assume  that  the  behaviour  of  the  dipole  density  at  bow 
and  stem  at  the  surface  is  the  same,  for  an  arbitrary  form,  as  for  a 
lenticular  cylinder  having  the  same  included  angle,  then  the  low- speed 
asymptotic  behaviour,  which  depends  only  on  certain  derivatives  evaluated 
at  the  bow  and  stem  at  the  surface  ([3],  [7]),  is  easily  derived  in 
either  theory,  allowing  comparison  of  the  theories  with  existing 
experimental  results. 


Incidentally,  it  can  be  seen  that  cannot  be  used  indis- 

hlr 

criminately  since  it  predicts  infinite  wave  resistance  for  a  vortical 


circular  cylinder.  The  theoretical  implications  of  this  divergence 
arc  being  studied.  The  possibilities  of  practical  application  of  both 
and  CTpjp  ^ire  presumably  limited  to  bodies  which  more  nearly  resemble 
a  ship  thaii  a  verticle  circular  cylinder  docs. 
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APPENDIX  A 

CONVERGENCE  OF  THE  SERIES  FOR  g^(x)  and  f^(x)  AS  FUNCTIONS  OF  € 
If  we  let  p  =  2-s,  then  from  Equation  (5) 


(A-^1) 


_  sir  -  sin  stt 
8  sin^ 


- t!  [1  +  •••) 


For  thin  cylinders,  when  2  ^  p  >  1,  we  have  0  ^  s  <1,  and 
0  ^  e  <  tt/  8 . 

Let  us  now  regard  s  as  a  complex  variable.  The  function 
CSC  ^  has  simple  poles  at  s  =  0,  d"  2,  +  A,  but  at  s  =  0  we  have 

c  =  0;  thus  €  is  an  analytic  function  of  s  for  Is  j  <  1,  the  region  we 
are  interested  in. 

Now  consider  g^(x) , 

(A-2)  g_(x)  “  Y  (  'X  cot^  ^  +  1  -  4x^  -  cot 


The  function  cot  has  simple  poles  at  s  =  0,  "t  2,  •••,  but 

'  2  ""  2 
cot^  ^  +  1  -  4x^ 


has  a  branch  point  when  s  =  s^,  where  cot' 


TTS 


o  =  -  (l-4x^),  or 


^  =  I  t  2n  +  -^  Arc  tan  ( i  \/  l-4x^  ),  n  =  0,  1,  2,  ••• 

Thus  1Sq|2  g^(x)  is  an  analytic  function  of  s  for 

si<  1 
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Thus  |s|^|  s  1,  and  is  an  analytic  function  of  s  for 

is'  <  I. 

Therefore,  both  g  (x)  and  f^(x)  can  be  expressed  as  a  convergent 

power  series  in  s  for  0  s  <1. 

For  s  0,  £  0,  s  and  £  are  linearly  related  (see  Equation  (A-1)), 

and  hence  if  a  convergent  power  series  expansion  exists  in  powers  of  s,  such 

an  expansion  also  exists  in  powers  of 
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APPENDIX  B 

RELATIONS  BETWEEN  NORMALIZED  AND  UNNORMALIZED  QUANTITIES 


Let  g  (x)  define  a  lenticular  cylinder  of  length  L  via 


=  + 


g,(x), 


AAA 

where  €L  is  average  half-beam.  If  we  let  z  =  zL,  y  ==  yL,  x  -  xL  and 
define  g^(x)  =  L  ^g^(xL)  we  find 


z  =  f  g^(x) 

as  the  normalized  form  of  the  lenticular  cylinder.  If 


^  .(x,z) 


A 

+  cx 


is  the  total  potential  describing  a  uniform  flow  past  a  lenticular 

A 

cylinder  in  unnormalized  coordinates,  and  if  we  let  ^^(xL,  zL)  =  ^g(x,z) 
we  find  that  the  total  normalized  potential  is 


^^(x,z)  +  cLx. 

The  l-Inearized  relation 

A  A  A 

d^^(x,o  ir)  dg  (x) 

— ^  sr- 


becomes 

dg^(x) 

dz  ■  ^  3x 

in  normalized  coordinates  Since  the  dipole  density  is 

A 

X  ^  A  1 

^  /  ^^^(x,0^)  A 

-  TT-  I  - -  dx  we  have  the  following  linearized  relations: 

^TT  J  ° 
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C  ^  ^  V 

Unnormalized  dipole  density  “  “ 

Normalized  dipole  density  “  “  ^  SgC^t)  . 

For  completeness  we  note  that  the  free  surface  condition  in 
normalized  coordinates  is 

^  =  F  ,  (y  -  0) 

dx  dy 

where  F  =  ^  . 

c 
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APPENDIX  C. 

SUPERCIRCULAR  CYLINDERS,  0  p  1 

A 

Consider  the  velocity  potential  i>  as  given  by  Equation  (4): 

A 

it  becomes  infinite  whenever  q  =  0  and  %  =  mpTr,  wliere  m  =•  1,  2,  3,  ••  •  . 
Since  in  the  region  of  interest,  z  ^  0,  we  must  have  0  £  ^  £  ir,  the 
restriction  on  m  is 
(C-1)  mp  1  . 

Tlicrcfore,  for  p>l  (subcircular  cylinders)  no  m  satisfies  Equation  (C-1). 
For  supcrcircular  cylinders,  p  =  p^  1,  there  are  m^  singularities  in  ^), 
vvhcrc  m^  is  the  largest  integer  satisfying  the  equation  m  ^/Po° 

The  rectangular  coordinates  corresponding  to  q  =  0  and 
^  -  mp-rr  arc 

X  =  0 

m 

(C-2) 

^  a  a  sin  mp-rr 
^m  1-cos  mp-rr 


Note  that  the  point  (0,z^)  for  m  ^  1  is  inside  the  cylinder  because 
i  =  mpm  >  • 

Lot  us  define  local  polar  coordinates  r,  a  with  origin  at 
(0,  z^)  ,  see  Figure  11.  As  q  ->  0  and  ^  ->  mpTr  Equation  (4)  can  be  v/ritten 
ns 

(C-3)  L  ->  2ca  -2 - 2 - 

™  r^^+imp-jT-O 

which  in  terms  of  r  and  a  for  r  «  1  becomes 
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(C-4) 


cosa 

r 


where 


(C-5) 


2c  a 

L-COS  mpTT 


Equation  (C-4)  is  the  potential  due  to  a  line  dipole  source. 

2 

For  p  =  1  (circular  cylinder)  m  =  1  and  =  ca  ,  “  0,  which  just 

expresses  the  well  known  fact  that  a  circular  cylinder  of  radius  a  can  be 

generated  in  a  uniform  flow  with  velocity  c  by  a  dipole  distribution  of 

2 

strength  =  ca  on  the  axis  of  the  cylinder. 

Thus  for  p  <  1  in  addition  to  any  dipole  density  -  fg(x) 
on  the  x-axis  we  also  have  discrete  dipole  sources  of  strength  at 
(0,  z  ).  As  p  decreases,  z  for  a  given  m  increases,  and  so  does  p  ; 
also  new  discrete  dipoles  are  created  at  the  origin  as  p  decreases 
(whenever  P  “  ~  some  m)  ,  which  then  split  and  move  away  from  the 
origin  in  the  negative  and  positive  z- direct ions  as  p  decreases  further. 

So  far  we  have  kept  the  length  2a  constant,  which  means  that 
the  circles  forming  the  cross  section  increase  as  p  0.  Now  let  us 
keep  the  radius  k  of  the  circles  constant,  in  which  case  2a  decreases 
as  p  0, 


(C'6)  a  =  k  sin 

Equations  (C-2)  and  (C-5)  become 

sin  sin  mpir 


(C-7) 


(C-8) 


A 

z  =*  k 
m 


L-cos  mpir 


=  9  ^ 

^  ^  l-cos  mpir 
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If  we  now  let  p  -♦  0  while  k  remains  constant,  the  circles 
become  tangent  at  the  origin,  and  there  is  an  infinite  number  of  discrete 
dipoles  on  the  z-axis  at 


(C-9) 


m 


=  12  3 
^  y  ^  y  y 


of  strength 

(C-10)  li,  =  % 

m 

Note  that  z^  never  exceeds  k,  and  the  strength  of  the  dipoles 
near  the  origin  approaches  zero. 
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FIGURE  lA 


LENTICULAR  CYLINDER  IN  A  SEMI-INFINITE  FLUID, 


UNNORMALIZED  COORDINATES 


LENTICULAR 


FIGURE  IB 

CYLINDER  IN  A  SEMI-INFINITE 
NORMALIZED  COORDINATES 


FLUID, 


SUBCIRCULAR  LENTICULAR  CYLINDER 


A 


FIG.  2B 
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